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EXPERIMENTAL  TECHNIQUE  AND  EXPERIMENTAL  METHOD 

PARALLEL  CRYSTAL  LASER  POSITION  SHIFT  INTEROKETER  AND  ITS 
APPLICATION  IN  FLOW  FIELD  MEASUREMENT 

Liu  Jianbang,  Xia  Shengjie  Institute  of  Mechanics, 

Academia  Sinica 

Abstract :  The  principle  of  the  instrument  is  to 

obtain  an  interference  pattern  by  reflecting  from 
the  front  and  rear  surfaces  of  a  parallel  crystal 
laser  beam  (after  beam  expansion)  after  passing 
it  through  the  system  to  be  measured.  The  instru¬ 
ment  has  good  vibration-preventive  property,  and 
can  be  used  conveniently  to  measure  flow  fields 
with  relatively  large  lateral  and  longitudinal 
dimensions.  The  instrument  has  high  light  source 
intensity  and  simple  structure,  and  may  be  easily 
calibrated  and  used.  It  has  a  good  prospect  of 
being  widely  used  in  flow  field  measurement, 
especially  in  the  measurement  of  supersonic  flow 
field. 


I .  Forward 

An  interferometer  is  an  effective  flow  field  measuring 
instrument.  The  commonly  used  Mach-Zehnder  interferometer 
(abbreviated  M-C  interferometer)  is  less  tolerant  of  vibra¬ 
tions,  has  smaller  field  of  view  and  depth  of  measurement. 

Its  light  source  has  lower  intensity  and  hence  requires  more 
stringent  demand  on  the  measuring  environment  with  more 
limitations  on  the  lateral  and  longitudinal  ranges  of  the  flow 
field  being  measured.  It  Is  also  difficult  to  make  continuous 
recording  of  high  speed  flow  processes.  The  reasons  for  these 
shortcomings  are:  1)  There  are  4  optical  elements  other  than 


1 


the  light  source  and  the  recording  element  that  are  involved 
In  the  interference  process.  This  creates  3  problems:  during 
exposure  there  cannot  be  any  relative  shifting  of  positions  of 
the  4  optical  elements  comparable  to  a  wave-length  (this  is  a 
major  shortcoming  of  interferometers  in  general).  The  method 
to  overcome  this  problem  is  to  use  vibration-free  base  which 
tends  to  make  the  instrument  bulky  and  heavy,  inconvenient 
both  to  manufacturing  and  to  utilization;  precise  calibration 
of  the  relative  positions  of  the  optical  elements  is  necessary. 
The  more  elements  there  are,  the  more  complicated  becomes  the 
calibration  process  and  the  more  stringent  is  the  demand  on 
the  design  and  construction  of  the  calibration  mechanism; 
optical  elements  require  precision  optical  work.  The  more  ele¬ 
ments  there  are,  the  larger  will  be  the  quantity  of  optical 
work.  This  is  even  more  so  for  making  an  interferometer  with  a 
large  field  of  view.  2)  Before  the  advent  of  the  laser,  light 
source  intensity  for  interferometers  is  low  with  inferior  qual¬ 
ity  of  interference,  not  suitable  for  the  need  of  high  speed 
photography  as  well  as  imposing  limitation  on  the  longitudinal 
range  of  the  measured  flow  field. 


In  this  paper  we  investigage  a  parallel  crystal  laser 

position  shift  interferometer  (abbreviated  as  parallel  crystal 

interferometer).  Its  principal  feature  is  that  only  one  single 

optical  element  is  used  to  produce  the  interferogram,  namely 

the  front  and  rear  surfaces  of  a  crystal  beside  the  light 

source  and  the  recording  element.  Another  feature  is  that  a 

laser  is  used  as  the  light  source.  Our  work  showed  that  the 

parallel  crystal  interferometer  has  overcome  in  a  great  measure 

M-C 

the  shortcomings  of  the ^  interferometer  mentioned  above,  and 
has  the  prospect  to  become  an  interferometer  with  high  utili¬ 
zation  value,  which  may  become  widely  used  in  flow  field 
measurement . 
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For  the  general  theory  and  application  of  phase  shift 
Interferometers  may  be  found  in  reference  [1].  A  parallel 
crystal  was  used  as  the  basic  element  of  a  phase  shift  inter¬ 
ferometer  in  references  [2]  and  [3]  but  the  method  in  [2] 
produced  interference  fringes  with  relatively  low  contrast 
while  reference  [3]  merely  quoted  the  results  of  [2].  In  this 
paper  we  analysed  the  principle  of  a  parallel  crystal  inter¬ 
ferometer,  explained  its  advantages  and  disadvantages  and  their 
causes,  quantitatively  analysed  the  formation  of  the  interfer¬ 
ence  fringes  and  cited  our  experimental  results. 

II.  Basic  Principles 
1.  General  discussion. 

As  shown  in  Figure  1,  a  laser  beam  with  good  direction¬ 
ality  and  coherence  is  expanded  into  a  plane  wave  with  a  larger 
cross-section  and  even  better  directionality  after  passing 

through  a  beam  expander.  After  travelling  through  the  system 

measured 

flow  parallel 

bean  expander  field  crystal 


observation  screen 

Figure  1.  Principle  of  a  parallel  crystal  interferometer 

being  measured,  it  is  duly  reflected  from  the^front  and  rear 
surfaces  of  a  parallel  crystal  to  form  two  beams  with  a  slight 
shift  in  position,  the  overlapping  portion  of  which  now  inter¬ 
fere  and  produce  interference  fringes  on  the  observation  screen. 
If  the  measured  field  flow  is  ideally  uniform,  then  an  ideal 


plane  wave  through  it  will  remain  an  ideal  plane  wave.  We 
shall  prove  below  that  an  ideal  plane  wave  will  produce  an 
interference  pattern  with  equi-distant  parallel  lines  after 
reflecting  from  the  front  and  rear  surfaces  of  a  parallel 
crystal.  Therefore,  from  the  deviation  of  the  interf erogram 
of  the  measured  flow  field  from  a  diagram  of  equi-distant 
parallel  lines,  we  can  make  direct  qualitative  judgment  on  the 
density  distribution  of  the  measured  flow  field  as  well  as 
obtain  quantitative  conclusions  from  computation. 

2.  Theoretical  analysis  of  the  interferogram 

We  assume  that  the  light  beam  propagates  along  the  z 
axis  in  a  Cartesian  coordinate  system  xyz.  (Figure  2). 

cl 

(x,y)  is  the  phase  distribution  function  in  the  cross-section 
area  A.  The  2  beams  reflected  from  the  parallel  crystal 
deviates  by  an  amount  Ax  along  x  axis.  Ax  is  called  the  dis¬ 
tance  of  position  shift,  and  x  axis  the  axis  of  position 
shift.  The  phase  distribution  functions  of  the  2  beams  on 
cross-section  B  are  respectively  ^Bp(x,y)  and  jZ<BR(x,y).  The 
subscripts  "F"  and  "R"  indicate  whether  the  reflection  is  from 
the  Front  or  the  Rear  surface.  After  reflection,  the  coordi¬ 
nate  system  propagates  along  the  reflected  beam  from  the  front 
surface.  Hence 

<£*«(*>  v)  =  <£.<(*,  y),  </>««(*> y)  =  4-  Ax, y)  +  a  +  bx  £2) 

where  a  +  bx  is  the  phase  change  due  to  reflection  from  the 
rear  surface.  It  is  not  difficult  to  find  the  constants  a  and 
b  from  the  crystal  structure,  its  position  and  properties  of 
its  material,  b's  value  obviously  is  directly  proportional  to 
the  angle  between  the  front  and  rear  surfaces  of  the  parallel 
crystal  (i.e.  its  wedge  angle).  On  cross-section  B,  (x,y) 

and  /BR(x,y)  has  a  phase  difference 


M 

cross-sue t ion.4  I 
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(flow  field 
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cross-section  B> 


Parallel  crystal 


Figure  2 


By  definition 
then 


^a(r>  y)  «  </>ac(*,  y)  —  </>b.s(.t,  y) 

“  +  Ax,  y)  ~f-  a  +  £x  —  y) 

A./>X*.  y)  =  4>A*  +  Ax,  y)  -  *,(*,  y) 


y)  -  A<t>Xx,  y)  +  a+  bx 


The  surface  with  equal  phase  for  an  ideal  uniform  flow  field 
is  the  xy  plane,  therefore  ti^(Xjy)  I  theory  =0  and  hence 

<>»<*,  y)  I  ■  •=  O  +  bx  (  4  ) 

A<#>X*,y)  —  ^B(x,y)  —  <£*(*,  y)  I  (5) 

Now  we  shall  find  the  phase  difference  4>a  (xK,  y)  —  <^(Xo,  y), 
of  2  points  (xo,y)  and  (x^,y)  in  cross-sectionj^A  where  the 
difference  between  x^  and  xQ  is  an  integral  multiple  of  the 
position  shift  distance,  namely  *n  —  x0  *=*  WAx,  n  —  i,  2,  •  • .. 

We  first  divide  the  line  segment  x^  -  xQ  into  N  equal  parts, 
the  coordinates  of  the  individual  points  of  division  being 
xn,  n  =  0,  1,  2,  N.  (Figure  3).  From  equations  (2)  and 

(5)  we  can  find 


4>a(*n>  y)  —  4>A*ot  y)  *“  2  y)  —  ^C***  y)l»l 
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Figure  3 

We  shall  now  discuss  the  above  results  as  follows:  1) 
equation  (*J)  shows  that  for  an  ideal  uniform  flow  field  the 
phase  difference  of  the  2  beams  reflected  from  the  front  and 
rear  surfaces  of  the  parallel  crystal  is  directly  proportional 
to  x.  Therefore,  the  interferogram  of  an  ideal  uniform  flow 
field  is  a  set  of  equi-distant  parallel  straight  lines  normal 
to  the  position  shift  axis,  the  distance  being  2ir/b .  From 
this  we  know  that  the  distance  between  the  parallel  lines  must 
be  inversely  proportional  to  the  wedge  angle  of  the  parallel 
crystal.  The  deviation  of  the  interferogram  of  the  realistic 
flow  field  from  that  of  the  ideal  uniform  flow  field  is  due  to 
the  non-uniformity  of  the  density  of  the  medium  in  the  corre¬ 
sponding  region.  From  this  we  can  determine  the  density  dis¬ 
tribution  of  the  measured  flow  field  without  too  much  diffi¬ 
culty  . 

2)  From  equation  (6)  we  can  quantitatively  determine  the 
difference  in  the  average  density  along  two  straight  lines 
parallel  to  the  optical  axis  in  the  measured  flow  field: 

p(xN,  y)  “  &(.*<»  y)  —  —  1  [#>(*w*y»0  —  p(*o*  y»  O  Vs 

I  J9 

The  projections  of  these  two  straight  lines  on  the  observation 
screen  have  coordinates  (x^,y)  and  (xQ,y)  while  the  two  end 
surfaces  of  the  measured  f]ow  field  are  assumed  to  be  the  z=0 
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and  z=l  planes.  The  process  is  as  follows:  Choose  the  coordi¬ 
nates  of  the  projections  of  the  observation  screen  of  the  2 
straight  lines  to  be  measured  to  be  (x  y)  and  (x  .y).  They 

should  satisfy  x  -x  =HAx,  M=i,2....  Rotate  the  crystal  about 
no 

its  central  axis  so  that  the  position  shift  axis  is  parallel 
to  the  line  connecting  the  2  points.  After  the  interferogram 
is  photograohed ,  read  off  from  it  the  value  of  2$n(X  .y)  for 
n=0, 1 , . . . , N-l ,  using  (xQ,y)  as  point  0.  Then  measure  the  dis¬ 
tance  between  the  interference  fringes  in  the  vicinity  of 
(xQ,y)  and  find  the  value  of  b  in  the  equation  (4).  With  (xQ,y) 
still  as  the  zero  point,  find  the  value  of  tfg(xn,y)  ideal  at 
n=0, 1> • • • j N-l  from  equation  (4),  then  from  equation  6  we  can 
find  the  phase  difference  between  the  points  (x^,y)  and  (xQ,y) 
on  cross-section  A.  Finally  find 


y)  —  p(x o>  y)  =  —rz  [4>a(,*n>  y)  —  <^(r°=  y)l 

2nlK 


(7) 


where  A  is  the  wave-length  and  K  is  the  Glaston-Derr  constant. (4) 


3)  Equation  (3)  shows  that  the  phase  difference  $g(x,y) 
of  2  light  rays  incident  on  one  point  on  the  observation  screen 
is  the  sum  of  the  phase  difference  Ad^(x,y)  of  the  corresponding 
2  points  on  the  output  cross-section  of  the  measured  flow 
field  (their  distance  is  equal  to  the  distance  of  position 
shift  Ax)  and  the  phase  difference  a  +  bx  of  the  2  light  rays 
after  being  reflected  from  the  front  and  rear  surfaces  of  the 

A 

parallel  crystal.  The  value  of,£iA(x,y)  for  various  measured 
flow  fields  along  the  direction  of  the  optical  axis  with 
different  average  densities  and  lengths  does  not  change  very 
much  because  the  light  ray  in  the  region  of  the  measured  flow 
field  has  not  yet  been  split  into  2  beams.  Magnitudewise , 
(A/2*)A</>.<(r,  y)  is  generally  of  the  order  of  the  wave  length, 
and  (1/2 jr)(«  +  £x)  is  generally  of  the  order  of  centimeters, 

i.e.  (a  +  $x)  »  y)  •  Hence  in  the  parallel  crystal  inter- 
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ferometer,  the  value  of  the  optical  path  difference  of  the  2 
interfering  beams  is  primarily  determined  by  (i/2i.)  («  t  ix) 
namely  by  the  characteristics  of  the  runnel  crystal  and  its 
value  is  of  the  order  of  several  centimeters.  Therefore  the 
parallel  crystal  interferometer  does  not  require  the  use  of 
the  compensator  that  is  commonly  used  in  M-C  interferometers, 
which  is  calibrated  according  to  the  length  of  the  optical  path 
of  the  measured  flov/  field.  Also  the  requirement  on  the  coher¬ 
ent  length  of  the  laser  light  source  is  not  very  stringent. 

3.  Advantages  and  disadvantages 

1.  The  principle  optical  elements  that  produce  the  inter¬ 
ference  effect  in  the  M-C  interf erometer  are  the  2  beam  splitt¬ 
ers  and  the  2  reflecting  mirrors.  Any  small  change  in  the  rela¬ 
tive  positions  of  these  4  elements  will  affect  the  optical  path 
difference  of  the  2  interfering  beams;  the  parallel  crystal 
interferometer  uses  only  1  optical  element  to  produce  interfer¬ 
ence.  Even  under  strong  vibrations,  the  stability  of  the 
relative  position  of  the  front  and  rear  surfaces  of  the  para¬ 
llel  crystal  remains  to  be  unquestionable .  This  is  the  reason 
why  the  parallel  interferometer  has  the  capability  to  with¬ 
stand  strong  vibrations. 

2.  The  measured  flow  field  is  placed  between  the  parallel 
beam  expander  and  the  crystal  where  the  beam  has  not  yet  been 
split,  therefore  the  relative  vibration  of  the  beam  splitter, 
the  measured  system  and  the  crystal  has  no  effect  on  the  opti¬ 
cal  path  length  difference  between  the  2  interfering  beams; 

as  shown  on  the  observation  screen,  the  field  of  view  of  the 
interferogram  may  have  considerable  vibrations  but  the  inter¬ 
ference  fringes  themselves  remain  stable.  In  practice,  the 
beam  splitter,  the  measured  system  and  the  parallel  crystal  may 
be  mounted  on  three  Independent  stands  which  require  no  special 
set  up  to  eliminate  vibrations.  This  greatly  simplifies  the 
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experimental  set  up  and  facilitates  the  measurement  of  flew 
field  with  relatively  large  longitudinal  (along  optical  path 
direction)  dimension  (it  is  entirely  feasible  to  measure  flow 
field  several  meters  in  length).  In  addition,  the  4  optical 
elements  that  produce  the  interference  effect  in  the  I'.-C 
interferometer  subtend  an  angle  of  45°  with  the  optical  axis 
so  that  their  diameters  should  be  times  that  of  the  light 

beam.  Two  compensator  windows  with  diameters  slightly  larger 
than  the  beam  diameter  are  also  required.  In  the  parallel 
crystal  interferometer  there  is  only  one  optical  element  to 
produce  the  interference  effect  with  its  angle  subtended  with 
the  incident  beam  nearly  perpendicular  and  hence  its  diameter 
needs  only  be  slightly  larger  than  that  of  the  beam.  There  is 
no  need  for  a  compensator  in  the  parallel  crystal  interfero¬ 
meter,  thereby  lowering  considerably  the  demand  for  optical 
work  and  moulding  for  optical  element  vis  a  vis  the  M-C  inter¬ 
ferometer  . 

3.  With  the  laser  as  the  light  source,  the  intensity  is 
much  higher,  making  it  possible  to  make  continuous  recording 
of  high  speed  processes  with  high  speed  cameras. 

4.  Simplicity  in  structure  and  convenience  in  calibra¬ 
tion  and  utilization  are  reflected  in  fewer  optical  elements, 
simple  supporting  stands,  relaxed  requirement  on  the  coherent 
length  of  the  laser  light  source,  the  elimination  of  compensa¬ 
tors,  the  lack  of  necessity  to  make  precise  calibrations  on 
the  relative  positions  of  different  optical  elements,  the 
directness  of  the  observed  results,  etc. 

5.  The  disadvantage  is  that  in  determining  the  average 
density  difference  of  2  straight  lines  in  the  flow  field 
quantitatively.  It  is  necessary  to  add  several  observed  results 
together,  making  the  procedure  more  complicated  than  that  for 
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the  M-C  interferometer;  for  any  non-transparent  object  in  the 
flow  field,  there  will  be  "double  shadow"  around  the  edge  of 
the  object  on  the  observation  screen.  Although  this  effect 
may  be  reduced  by  shortening  the  position  shift,  yet  in  gen¬ 
eral  it  cannot  be  completely  eliminated.  We  must  point  cut 
that  there  is  no  double  shadow  in  the  interference  fringes  so 
that  the  effect  of  the  double  shadow  generally  will  not  affect 
the  result  of  measurement.  It  will  only  give  some  uncomforta¬ 
ble  feeling  to  the  observation. 

III.  Experimental  set  up  and  experimental  results 

The  parallel  crystal  interferometer  may  be  subdivided  into 
3  parts:  the  parallel  light  source,  the  parallel  crystal  and 
the  recording  device.  The  parallel  light  source  consists  of 
the  laser  and  the  parallel  beam  expander.  An  ordinary  single 
mode  continuous  output  or  pulsed  output  laser  tube  may  be  used 
for  the  laser.  For  the  parallel  beam  expander,  we  may  use 
either  the  transmitting  or  the  reflecting  type.  Its  effect  is 
to  expand  the  beam  diameter  to  the  desired  size  and  to  reduce 
the  angle  of  divergence  to  obtain  a  better  plane  wave.  The 
parallel  beam  expander  should  be  focused  to  infinity.  The 
diameter  of  the  beam  incident  on  the  object  lens  should  be 
slightly  larger  than  the  diameter  of  the  lens  so  that  only  the 
central  portion  of  the  beam  is  transmitted  in  order  to  warrant 
a  more  uniform  intensity.  The  position  of  the  laser  and  the 
beam  expander  should  be  suitably  arranged  so  that  the  hot  air 
flow  produced  by  the  laser  tube  will  not  affect  the  quality 
of  the  beam.  The  laser  and  the  beam  expander  should  be  mounted 
on  the  same  base  support  to  reduce  the  movement  of  the  field 
of  view  induced  by  vibration.  We  used  as  light  source  a  single 
mode  continuous  output  HeNe  laser  with  a  power  rating  of  1  mW. 

A  reflecting  optical  path  beam  expander  was  used  with  an  eye¬ 
piece  of  focal  length  15  mm,  and  an  objective  of  focal  length 
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1  m  and  diameter  120  mm.  The  output  beam  was  co°cl  with  uni¬ 
form  intensity.  The  structure  of  the  beam  expander  and  its 
support  is  shown  in  photo  1  on  picture  page  1. 

The  front  and  rear  surfaces  of  the  parallel  crystal  should 
be  good  quality  optical  surfaces  without  any  coating.  The 
wedge  angle  should  be  between  a  few  seconds  and  a  few  ten's  of 
seconds.  There  is  no  strict  requirement  except  to  produce  an 
interference  pattern  with  suitable  distance  between  the  fringes. 
We  used  a  wedge  angle  of  6".  The  first  and  second  reflected 
light  from  the  front  and  rear  surfaces  of  the  un-coated  crystal 
had  nearly  equal  intensity  while  that  of  the  third  and  higher 
reflection  are  rather  weak,  producing  an  interferogram  with 
good  contrast.  There  was  practically  no  effect  on  the  inter¬ 
ferogram  from  the  3rd  and  higher  reflected  light.  For  example, 
if  the  parallel  crystal  is  made  of  K^glass  with  an  index  of 
refraction  of  0.042,  then  the  ratio  of  the  incident  light 
intensity  to  the  first  few  reflected  and  transmitted  light 
intensity  are  shown  in  Figure  4.  If  is  easy  to  show  that  in 
the  interferogram  of  the  first  and  second  reflected  light  the 


Figure  4.  Intensity  ratio  of  incident  light  and  the  first  few 
reflected  and  transmitted  light  for  a  glass  crystal. 

intensity  of  the  darkest  spot  is  only  0.006  times  that  of  the 
brightest  spot.  Hence  the  bright-dark  contrast  of  the  inter¬ 
ferogram  is  close  to  ideal.  The  transmitted  part  is  about  92!?, 
which  may  be  used  for  other  measurements. 
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The  recording  may  be  done  by  direct  observation  on  the 
screen  or  by  still  camera  or  high  speed  movie  camera.  The 
observation  screen  is  a  flat  white  diffuse  reflecting  surface. 

A  piece  of  ground  glass  may  also  be  used.  When  the  fringes  are 
too  crowded,  an  enlarging  device  may  be  used  in  front  of  the 
screen.  If  a  still  camera  or  a  high  speed  movie  camera  is  used, 
it  may  be  mounted  without  the  lens  in  place  of  the  screen  to 
record  the  interferogram  directly  on  the  film. 

The  three  parts  mentioned  above  were  mounted  separately  on 
3  independent  stands  (photo  1  picture,  page  3).  The  stands 
have  no  special  vibration-isolating  device. 

After  the  proper  installation  of  the  apparatus,  an  inter¬ 
ferogram  will  appear  on  the  screen.  Since  now  the  observed 
region  contains  uniform  density  air,  the  interferogram  should 
consist  of  equi-distant  parallel  straight  lines  perpendicular 
to  the  position  shift  axis.  From  this  we  can  examine  the  qual¬ 
ity  of  the  beam  expander  and  of  the  crystal.  According  to 
experimental  need,  we  may  rotate  the  crystal  about  the  central 
axis  to  obtain  the  required  fringe  direction  and  then  about  the 
y  axis  to  calibrated  the  distance  of  position  shift.  Finally 
the  recording  device  is  to  be  calibrated.  Except  the  recording 
device,  the  whole  calibration  procedure  may  be  accomplished  in 
a  few  minutes. 

The  experimental  results  are  as  follows:  Photo  2,  picture 
page  5  shows  the  interferogram  of  uniform  density  air  in  the 
region  of  measurement.  It  approaches  very  closely  the  equi¬ 
distant  parallel  straight  lines  normal  to  the  axis  of  position 
shift  as  calculated  from  equation  (*0.  Photo  3,  picture  page  6 
is  one  of  the  interferograms  oblique  laser  beam  recorded  con¬ 
tinuously  with  a  high  speed  movie  camera  for  a  super-sonic 
experiment  under  strong  vibration.  The  precision  of  this  type 
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of  measurement  depends  on  the  method  to  measure  the  distance 
between  the  fringes.  Photo  ^  in  picture  page  11  is  the  inter- 
ferogram  of  the  warm  air  flow  field  near  a  human  hand.  From 
the  manifest  curvature  of  the  fringes  at  the  palm  we  can  see 
that  the  parallel  crystal  interferometer  is  relatively  sensi¬ 
tive  . 

Because  of  the  position  shift  of  the  two  reflected  beams 
due  to  the  thickness  of  the  crystal  and  the  incident  angle, 
double  shadow  appears  in  the  interferogram .  The  thicker  is 
the  crystal  and  the  larger  the  angle  of  incidence,  the  more 
separation  is  there  in  the  double  shadow.  To  overcome  this 
disadvantage,  we  shall  reduce  the  thickness  of  the  crystal  so 
that  the  two  shadows  will  approach  each  other  and  overlap. 

For  this  purpose,  we  adopted  an  interference  element  (Figure  5) 
consisting  of  a  prism  and  a  plane  mirrow.  The  prism  is  in 
contact  with  the  mirror  with  a  thin  layer  of  air  between  them. 

The  parallel  beam  is  incident  on  the  prism  after  traversing 
the  measured  system  at  almost  right  angle  and  is  reflected  at 
the  two  surfaces  in  front  and  at  rear  of  the  air  gap  and  then 
emerges  from  the  other  surface  of  the  prism.  Interference  is 
produced  due  to  the  phase  difference  of  the  light  reflected  by 
the  two  surfaces.  This  basically  is  a  multiple  interference 
and  is  different  from  the  parallel  crystal  interferometer  dis¬ 
cussed  above.  Since  the  air  gap  is  very  thin,  the  position 
shift  of  the  reflected  beams  from  the  two  surfaces  is  also  very 
small.  If  the  air  gap  thickness  is  of  the  order  of  a  micro¬ 
meter,  then  the  position  shift  is  also  of  the  order  of  a  micro¬ 
meter  and  the  fuzziness  of  the  profile  of  the  model  in  the 
interferometer  will  also  be  of  the  order  of  a  micro-meter. 
According  to  photographic  theory,  a  picture  is  sharp  when  the 
diameter  of  the  fuzzy  element  is  less  than  0.1  —  0.05  mm. 

The  double  shadow  in  the  profile  of  the  profile  of  a  model  cannot 
be  discerned  in  the  interferogram  when  the  position  shift  is  of 


the  order  of  a  micro-meter.  Thus  the  phenomenon  of  the  double 
shadow  has  been  eliminated.  At  the  same  time,  fringes  with 
better  contrast  and  finer  details  than  that  found  in  the 
interferograms  of  apparallel  crystal  interometer  are  produced. 
Furthermore,  the  reflection  at  the  bottom  surface  of  the  prism 
is  almost  a  total  internal  reflection,  the  brightness  of  the 
interf erogram  is  much  greater  than  that  of  the  parallel  crystal 
interferometer.  This  in  particular  is  of  advantage  to  the  high 
speed  photographic  recording  of  fast  processes  with  a  much 
shorter  exposure  time.  From  photo  5  and  6  in  picture  page  1^ 
we  can  make  the  comparison.  The  former  is  the  photograph  of 
a  hot  nail  taken  with  a  parallel  crystal  interferometer  while 
the  latter  is  the  photograph  of  the  same  phenomenon  taken  with 
the  improved  apparatus.  We  can  see  that  the  double  shadow  is 
eliminated,  and  that  the  brightness  has  increased  and  the 
fringes  finer.  This  is  suitable  in  the  measurement  of  such 
flow  field  as  that  of  a  supersonic  wind  tunnel  which  contains 
models . 


measured  flow  field 
parallel  light  beam  i  ! 

i  i _ 

*  i 
I  I 

L - J 


Figure  5*  Air  gap  interference  element 


IV.  Conclusions 

The  parallel  crystal  interferometer  is  better  suited 
for  measuring  many  types  of  flow  fields,  such  as  the  flow 
fields  in  aerodynamics,  fluid  dynamics,  physics  of  explosion, 
heat  conduction,  and  laser  physics;  it  can  be  adapted  to 
adverse  measuring  environment  and  has  greater  practicality.  It 
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is  simpler  to  design  and  construct  and  is  easier  to  be  call 
brated  and  used.  It  is  specially  good  for  measuring  super¬ 
sonic  flow  field.  It  is  also  good  for  testing  the  quality 
optical  elements. 
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RESEARCH  MOTES 


NUMERICAL  SOLUTION  OF  NAVIER-STOKES  EQUATIONS 

Guo  Benyu  Shanghai  University  of  Science  and  Technology 
Guo  Benqi  Jiujiang  Refinery 

I.  There  are  many  methods  to  compute  Navier-Stokes 
equations  but  few  are  able  to  prove  their  estimation  of  errors, 
mainly  because  it  is  difficult  to  treat  the  non-linear  term 
and  the  ratio  of  pressure  to  density.  Artificial  state  equa¬ 
tion  was  used  in  references  [3]  and  [4]  while  Guo  Benyu  [6] 
combined  the  method  in  reference  [5]  and  proved  the  stability 
of  the  implicit  scheme.  However  he  has  not  proved  the  conver¬ 
gence  of  the  explicit  scheme  nor  has  he  been  able  to  derive  the 
convergence  of  P  from  this.  Another  method  to  solve  the  Navier- 
Stokes  equations  is  to  compute  directly  the  Poison  equation  of 
P.  For  convenience,  let  us  take  the  following  two  dimensional 
equations  as  an  example: 


where  U,  V  are  respectively  the  velocity  components  along  x, 
y  directions  is  the  coefficient  of  viscosity,  f^(i-l,2,3) 
are  known  functions.  Furthermore  the  initial  and  boundary 
values  of  U,  V  and  P  are  assumed  known. 


In  this  paper  we  extend  the  result  in  [5  ]  and  construct 
several  schemes  which  are  convenient  for  computation.  We  also 
have  extended  a  non-linear  inequality  in  [  5  ],  thus  solving 
the  problem  of  rigorous  error  estimation. 


II.  Let  h,  t  respectively  denote  the  step  lengths  alone 
the  space  and  time  directions  of  the  net.  i  —  rjh1  ,  Q  is  a 
net  point.  R^  represents  the  collection  of  internal  points  of 
the  net  region  to  be  solved  for  with  as  its  boundary.  w(k) 
or  w  both  denote  w(ih,jh,kT).  -  „  «.  +  1), 

all  have  the  same  meaning  as  in  reference  [5].  We  further 
define 

(tv,  «)  “  S  ^W£>)«(0)»  II “'ll1  “  (“'j  Ill’ll  *=*  ll«'xl|j  + 

/(«0  *“  2 

.  P«*i+R, 

where  R^,  R^  are  the  collections  of  the  internal  points  with 
minimum  x,y  coordinate  values  respectively. 

/,(«/)  -  [Mlf  +  /(«/) 

Lemma  Let  w(k)  be  a  non-negative  net  function,  gQ,  L,  t,  d 
are  positive  constants.  p(T)  is  a  non-negative  continuous 
function  of  T,  and  when  Wj  <  d(i  --  1)  ,  the  functionals 

Hu.( k  —  O  “  ,  w(k  —  01  <  Ltv(K  ~  O 

c.(0  -  GU-(O),--*,  «'C01  > 


then  if 


and  when 


g»tv(Q)  ^p(T)  <  e~LT/‘‘  ■  dgs 


then 


Note  1. 


G„(0  <  p(  r>  +  r  2 


w(0  <  Bill  eLT/i, 


HM  -  C.  [,«>  +  =3ij]  +  [c,  f  c,  -  c,]  KO 

cw(0  “  U  —  c,t  —  u/(^  —  ovu> 

where  v(k)  is  a  non-negative  net  function  and  t,  h  are  suffi¬ 
ciently  small,  and  there  exists  constant  M  such  that  C^M  - 
<  C^,  then  we  may  take 

d  -  Mh\  gt  -  j,  L  -  CoCM  +  1) 
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Thus  when 


e-*V*"+o  •  d  Bj , 

2 

«/( K)  <  2p,*:»T(w+» 


Note  2.  In  Note  1,  p^Mth}T,  C, ct*=b}r  where  b^,  b^* 

are  constants,  r  is  to  be  determined,  Co(r)  is  a  continuous 

function  of  r, hence  if  T  <  T  for  T  sufficiently  small  so  that 

o  o 

2  T(Af0e,c«T»  <  1  ,  then  we  may  take 


Thence 


-r~h -»  ^  =  2C„  ft-i. 

b\  •¥  bi  2 


»(ij)  ^  2pe*c°7 


III.  The  double  scheme  to  solve  ( UL )  is 

«,  +  */(«> *0  —  vA(«  4-  <rr«i)  4-  P£  **  /»>  Rh  X  S% 

«/,  4-  »,  *0  —  4*  art*,)  4-  P$  ”  /a>  -R*  X  5^  ► 

AP  ■»  «/>(«,  «0  4-/3,  R*  X  Sj 


(3.D 


where  0<<r,  a<l  are  non-negative  parameters,  denotes  the 
set  E{/|i  —  {tA  — 


If 


then 


<*>(«,  «0  -  («;*)»  “  +  ~ 

</(«/,  «,  «0  —  <*(«•«);  +  aO*0»  +  (1  —  a>«'i  +  (1  —  *)vu> 

a  -  1/2, «lr*  =  Hr*  -  0 

(«,  </(«,  u, «/))  -  (*,  “>  *'))  —  o 


This,  in  practice,  simulates  the  following  conservation  law  and 
hence  yields  a  better  numerical  result: 


iK"  f  + "  f ) t""'1  -  Sj(°  ^ 


V  Vdxdy  -  0 
dy  / 


Let  us  denote  the  errors  by  >  and  iet  e  — 

-  flsif*  +  BfC*. 
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then 


B*.D,“ll«ill,+  H?#0,»  ll£(OHoe.'=’l|f(Oli,+TV  2  Ii«CO!ii*  ll«,!!i=lli,!!i+l!?,!i; 

> -0 

Theorem  I.  In  scheme  (3-1),  if  (1)  •*  >  0,  (2>  > 

or  *  <  9v(i/2  —  a)’  *  ii'ere  exists  positive  number  Tq  and 

C  such  that  when  T  <  T  and 

o 

p  -  II  £  (  0  )  1]  *  +  rJ)S  i7.(;)l]'  <  V  •  C 

tm  0  i"1 

IU(0!I?®»  ^  Wp'tT»  where  L,  M  are  positive  constants  related  only 

to  m_,  and  u,  v  and  the  upper  limit  N  of  the  absolute 

v  2 

value  of  their  first  order  difference  quotients.  If  p  =  0(h  ), 

then  Tq  is  arbitrary. 

Proof:  Obviously  we  have  the  error  equations 


£,,(£,  5»  tt,  +  rf(2,  K  +  tt.  f +l')+li(#(  M,  *0  1'A(m  +  CTTK,)  *“  7l 

Lt(v ,  ii,  p")  “  £,  +  ^(2,  m~tm,  tf  +  tO  +  ^O'*  «>  *0  j/AC? -rffr?,) +  p>  =  7j 

Ap  “  </>(»,  ?)+?>  =  <#>(»>  «0  +  <£>(“»  SO  *t  *0  +  ?> 


(3.2) 


If  we  compute  the  quantity  (2m  +  mrS, ,  1,(2,  2,  £))  -§-  (22  +  mr2;,  ba(«,  ?,  ?))» 
then  similar  to  (5)  we  get 

Hell?  +  (m  —  l)rHe,(lJ  +  w»<e)  +  v  (racr  ~  \  ~  <T)S'(£‘'> 

+  tv  -r  *i(e)i  +  (v  —  «.*)IUU  —  ^  R(f) 

where  A,  B,  are  positive  numbers  unrelated  to  h  and  r  and 
a,  a^  are  positive  parameters  to  be  suitable  determined. 

8(e)  +  ■^-)c!le||,+l|e||,)+^-(IUI!,+ll6ll,ll£ilO+l!7ili,+ll7i!!:}+«,llp!l? 


Following  the  method  in  (5)>  we  choose  a  and  m  suitable  and  sum 
the  above  equation  in  the  t  direction  and  obtain 


IdOfli.'  <  M.B«C0)P  +  (*  +  j)"U<0)1  +  *•  2  *(e) 


(3.3) 
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From  the  third  equation  of  (3.2),  we  get 


where 


but 


'i(f)  "  (p>  J’)  +  (?,  7j)>  — t  2man7*!I*  +  l(p>/') 

hm0 


I  «•!  P 

»».=•  *«P  £-5 


|(?,  </■)!  <  ®  +  48  {MM?  +  iVJ||e||(  +  W,i!e:i»} 


thence 


||p!|{  <  96  {MM  +  a»!|£;ij  +  Ap||«|p]  +  4m$t* 


Substitute  the  above  into  (3.3),  then 


TlUCOliJ®.  -  -  AfJi||e(^)ll,l!e(^  -  Oil2  <  +  Mjr§  CD«C#>U* 

+  lleOW1)  +  r  53  W«(«i  MmMfU  +  0tQeO>n$  "^ileCOl!?) 

1*0  '  h 

Now  we  need  only  suitably  choose  and  apply  note  2  in  the  lemm 
in  which 

tv  “  “  II® lli»  r  “ 


Thus  if  U,  V,  P  are  suitably  smooth  so  that  (3-1)  has  an 
aymptotic  accuracy  of  at  least  0(h),  then  scheme  (3.1)  conver¬ 


ges  and 


OSiP  +  II? iP  =  o(a4),  ll?!P  =  o(a:) 


IV.  The  triple  scheme  to  solve  (1.1)  may  be  constructed 
as  follows: 


u,  +  0r2u,i  +  </(«,  u,  v)  —  vA u  +  pi  =  f,,  R»  X  5? 

v,  +  0t2v,\  •+•  d(.v,  u,  —  vAt>  4-  pt  “  /j,  J?*»X  S? 

A p  —  ^>(«,  »»)  4-  /j,  JJ*  X  S? 


(4.1) 


in  which  6  is  a  positive  parameter. 


k 


Theorem  II.  In  scheme  (*1.1)  if  0)»'>o,  ,  (2)  t9  is  constant. 


0  <  t6  < 


— ,  (3)7.  <  —  then  there  exist  positive  constants  T 


C  such  that  when  T  T„  p  —  p  -f  ||f(i)j;»  <  ch1 

IU(0!ifc*  +  l!«(*-l)i;,<  MpeL1 

2 

If  p'  =  0(h  ),  then  Tq  is  arbitrary. 


Proof:  Let  d  =*  ~  +  ii 

2r  t 

to  Theorem  I,  we  have 


,  then  -l/2<e„<0,  .  Similar 


+  (,»  -  +  ».,<«)  -  +  <»  -  -OB«« 

-  .  B IUD!  +r$»!!eB;i— 0»rl|ei!!,“0.rl!ei!!,+  H*<Bi  +  <  ^Ce) 


where  x(£)  is  the  same  as  in  Theorem  I.  Now  we  take  m  =  2/ 
2+*16o~9Av,  then  after  summing  over  t,  the  sums  of  the  2nd, 
4th  and  I2th  terms  in  the  LHS  of  the  above  equation  are  non¬ 
negative.  The  rest  of  the  proof  is  similar  to  that  used  in 
Theorem  I  but  we  take  w(0  “  ffeCMii©,  +  [j£(^  —  i)|j>  when 

applying  lemma.  ’ 


The  ADE  scheme  to  solve  (1,1)  is 


U,  +  d(^u,  u,  v )  —  vAh  —  —  W(«)  +  p:  *»  f>>  R»  x  5? 
v,  +  2(r,  «,  «')  -  vkv  —  ~  H(r)  -r -  Pi  =  1t>  R>  x  ^ 


A/>  —  <#>(«,  v)  +  /„ 


Rjxs;1 


(4.2) 


where  H(w)  assumes  the  forms  —  n>f,  —  uj,,  w„  —  «?i»  —  w'i.  +  **'«  + *•> 

etc  in  accordance  with  a  definite  rule. 


Theorem  III.  In  (scheme  (4.2)  v>0,  l<— ,  then  Theorem  I  is 

also  true. 


Proof:  Following  Theorem  I  we  have 

Ben;  +  On  -  i)r!!£,|!,+  —  *,0),  +  »•'.(£)  +  (v  -  a'0!!£H«  - 


-  ^  (£,WU)>  -  /(£,)  <  *(*) 

A  2 
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where  R(e)  is  same  as  before.  The  sun  of  the  last  two  terms  on 
the  LHS  is  not  more  than  ~  ||fi!{-f-2i’r  (6  +  ^)  !I£*U*  •  Now  if  we 

take  «  —  —  mv  then  the  sum  of  the  2nd,  7th  and  8th  terms 

2  —  i.v 

on  the  RH5  of  the  above  equation  is  non-negative.  The  rest  of 
the  proof  is  the  same  as  in  Theorem  I. 


The  splitting  scheme  to  solve  (1.0  may  be  constructed  as  fellows 


- - -  4-  — -  («*  +  u,  u) - — («?f  +  "it  J  +  “*  ft/^i  R*  X  $ 

r  2  2  2 


Lzift!  +  i2L  (s  +  „)  -  ^  is,5  +  «;?]  +  ?  -  nn,  R*  *  sj 

t  2  2  2 


f  4-  1^5-  (f*  +  v,  «)  —  — ■  [v*i  +  fiil  +  ^  =  /i/2,  R*  X  S? 


Ln^!  4-  *3-  c»  +  «0  -  ~  l?,i  +  M  +  ?  -  />/2>  R*  x  « 

r  2  2  ^ 


A/>  -*  <K"»  (;)  +  /<> 


Rt  X  S°r 


where 

dx(sw9  h)  *=  +  0 

</>•(«.,  v )  *=*  a(/^r)$  +  (1  — 

If 

a  —  “Mr*  “*  0 

2 

then 

(m-,  dx(.w,  «))  =  (w,  </)'(«>, 

(*1.3) 


Hence  the  conservation  law  is  simulated.  If  we  only  need  to 
compute  the  P  value  at  some  moment,  then  the  method  to  cal¬ 
culate  vorticity  and  flow  function  in  reference  [5]  may  be 
used.  Then  P  may  be  found  from  the  following  equation 


ylP 


+  /. 
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Similar  results  for  3  dimensional  problems  may  be  established. 
Further  results  may  be  seen  in  (6)  and  (7).  In  this  paper 
the  generalised  stability  index  of  the  corresponding  2  dimen¬ 
sional  scheme  in  reference  [6]  is  improved  to  s  <  1.  If  the 

reverse  wind  method  in  (6)  is  used  to  compute  (1.1),  the  scheme 
is  also  convergent. 
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